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We show that, contrary to common behef, internal shocks can arise in an accelerating 
radiation-dominated jet if it is confined even weakly to a converging opening angle be- 
cause the acceleration declines. The radiation-dominated internal shock (RDIS) enables a 
very efficient yet highly nonthermal emission by Fermi-like photon acceleration, keeping the 
electron-positron (e*) pair photosphere and inertia up to a high Lorentz factor > 1000. In 
gamma-ray bursts (GRBs), a weak confinement would persist beyond the progenitor star or 
surrounding matter because of the fast cocoon accompanying the breakout jet. The sim- 
plest model predicts few high-energy cosmic rays and neutrinos, and a correlation between 
the early afterglow and the GeV-TeV prompt emission. The central engine allows a less 
fine-tuned baryon load than previously thought, even including pure-leptonic unmagnetized 
outflows. 

\GRB and baryon load - Gamma- Ray Bursts (GRBs) are the most luminous objects 
in the universe. It remains a big challenge to reveal how most of the energy can be 
converted into gamma rays with highly nonthermal spectra.'^''^ 

The baryon load is a key parameter for the fireball dynamics and emission. 
The original fireball modeP'® was pure leptonic with the photon-to-baryon ratio 
r] = L/Mc^ = oo. As the fireball expands, the pairs annihilate and almost 
all the energy is released to thermal radiation: it is efficient but inconsistent with 
nonthermal observations. This leads to the baryon-loaded model with rj < 10'^, in 
which the radiation energy is converted to the kinetic energy of matter,^!' and then 
back to the radiation by internal shocks within variable outflows.'^^Sll gy its very 
nature, however, this standard picture faces an efficiency problem.l^''^^ The Lorentz 
factor dispersion required for the efficient reconversion is too large to reproduce the 
observed spectral correlations Vpeak oc ^1/2 xhis brings forth the photosphere 

modeP^"^ that invokes a partial thermalization near the photosphere to stabilize 
the peak energy Vpeak at the thermal peak, while retaining the nonthermal spectrum. 
The observed Band spectrum could be reproduced by Comptonization of the thermal 
photons via dissipation such as shock^lSli^2D an^j magnetic reconnections.'^''23 

These matter-dominated models assume the baryon load in the range (see Fig. [T]) 

to avoid the compactness problem^-* (for 10^ < rj) and the predominant thermal 
emissiorpS' (for rj < 10^), where the baryon-poor fireball with rj > 10'^ is considered 
to emit thermal photons because no internal shocks occur in the radiation-dominated 
fireball, which accelerates as -T oc r with the preceding shell always faster than the 
succeeding one. Therefore, the inferred baryon load is only limited within one order 
of magnitude, in contrast to a huge variety of GRBs. We may identify it as the 
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Fig. 1. We categorize GRB models on the basis of the baryon load. Popular models usually assume 
the photon-to- baryon ratio rj = L/Mc? from 10^ < r] (to avoid the compactness problem) to 
rj < 10^. We consider the radiation-dominated jet with r) > 10'', including the pure-leptonic 



original fireball model. 
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Such jets were considered to cause no internal shocks before being 



transparent at the Lorentz factor F ~ 10 and radius r/ro ~ 10 , and hence, emit (dominant) 
thermal radiation unlike in observations.'^ This picture is changed by the jet confinement. 

^^fine-tuning problem of the baryon load " . Although the baryon-rich fireballs could 
be just unobservable with selection bias, there seems no compelling reason to believe 
r] < 10^. Recent Fermi observationd^lUlS) have made the problem even worse^-^' by 
detecting high-energy photons to limit the Lorentz factor F > 10^, i.e., rj > 10^ in 
some bursts (despite the controversies on the numerical factoil^^^). 

One possible solution to the fine-tuned baryon load problem is the Poynting- 
dominated model.'^''^''^'^ Even without (electrons associated with) baryon, the 
energy is released only nonthermally via magnetic reconnections. Many uncertainties 
still remain, especially in the dissipation site. It may also be difficult to keep the jet 
cold in the hot photon pool near the central engine. 

The other possibility may be the externally induced dissipation by the oblique 
shock^21)-[36li Qj, ^Yie entrainment of baryon^^' during the acceleration. However, the 
external dissipation only occurs at the causally connected periphery of the jet, and 
the strong confinement {9 < r~^) is necessary for the full dissipation across the jet. 

In this paper, we show that, contrary to popular belief, the radiation-dominated 
jet can cause internal shocks if the geometry is confined even weakly (even 6 > r~^) 
by some external pressure or internal magnetic field. This could markedly expand 
the available parameter space to 10^ ^ as outlined below. 



§ Two shell model - We envisage that the Lorentz factor F of an outflow fluctuates 
by a factor of AF ~ a few around its mean. This may be modeled by two shells 
ejected from a radius tq at times ti and t2 = ti + At (where At > 0) with different 
initial Lorentz factors Fi and /2(> Fi), respectively (allowing for the possibility of 
Fi ^ 1 as in the jets breaking out of the progenitor star). The shell trajectory is 
given by 
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For the conical (or spherical) radiation-dominated jet F (x r, we can easily integrate 
Eq. © as r = {{r^/r^) + [c{t - tt) + roA]'}V2, (f = 1,2). 
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if the following condition is satisfied, 
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Two shells collide at 



(3) 



(4) 



If the interval At exceeds a critical time, cAt > (/32 — f3i)ro, i.e., 7^ < 0, even a 
photon cannot catch up with the accelerating fore shell (see Fig. [2]). 

The catch-up condition (jl]) is not usually fulfilled by the physical processes 
because the radial velocity dispersion spreads the shell width before the start r < tq 
and the next shell is delayed by the causal time. At > ro/c/3i for the non-relativistic 
case and At > rQ/2cr^ [and, hence, TZ < — /^^/(rf — F^)] for the relativistic case, 
leading to no collision. This is the essential reason why we have not considered the 
internal shock in the radiation-dominated jet, although a short timescale may be 
produced by a small-scale magnetic reconnection or a jet opening angle of less than 

~ i/r\ 

However, the Lorentz factor evolves differently from the conical one F oc r to 
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if the jet is confined to a converging opening angle, 

9(xr^-\ (A<1) 



(5) 



(6) 



by some external pressure or internal magnetic field. This is because the comoving 
volume evolves differently from the conical on^ to V ~ TrO'^r^FA oc r'^^F, where 
the lab-frame width A remains constant if it is smaller than the spreading width 
~ r/F'^. The adiabatic expansion of radiation (with an adiabatic index 7^ = 4/3) 
decreases the comoving temperature as T' cx V^~'^'^ oc V'~^/^ oc r' 

-2A/3^~l/3, A solid 

boundary conserves the total energy FT'^V' oc FT' = const, which yields Eq. ([5]). 

For the confined jet A < 1, the catch- up condition differs from Eq. We can 
analytically integrate the trajectory in Eq. ([2]) considering (hereafter) the relativistic 
case F2 > Fi ^ 1 with an expansion /3 ~ 1 - F~^/2. Then, the shehs colhde at 
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if Tlx > 0. 



and dZD. 



(7) 



Note that 7^A = 7^ for A = 1 and r2 > A » 1 in Eqs. 
The qualitative properties change at A = 1/2. 

(i) Strong confinement case, A < 1/2: The collision always occurs {TZx > 0) for 
any interval At, simply because the causal region ~ r/2F'^ oc r^^^'^ monotonically 
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Fig. 2. A spacetime diagram for the retarded time t ~ r/c in units of the causal time ro/2cri as 
a function of the radius log(r/ro), using F ^ 1. A photon travels horizontally in this diagram. 
A conical radiation-dominated jet F oc r (solid line) approaches a light path, so that any rapid 
jet emitted after the causal time ~ ro/2cFf cannot catch up with it. Meanwhile, a confined jet 
accelerates more slowly _r oc (thick solid line) in a converging opening angle 6 oc r^~^ with 
A < 1, enabling internal shocks with rapid jets (arrows) over a wide radius range. 



increases in the lab frame. This is almost the same as the usual internal shock 
after the acceleration ceases.'^® Meanwhile, since 9 < r~^, the jet boundary is 
also causally connected with the jet axis. This enables the recollimation or oblique 
shocks across the jet. Therefore, the jet can be dissipative during the expansion with 
shocks in the causal region. Note that the strong confinement {OF < 1) is employed 
to accelerate the Poynting-dominated jet to the matter-dominated jet.fS^'EB 

(ii) Weak confinement case, 1/2 < A < 1: The catch-up condition TZx > in 
Eq. d?]) requires that the interval At is less than a critical time Z\tcrit) 

This can be satisfied by the causal processes (i.e., with the next shell ejected after 
the causal time At > ro/2cr'^^), in contrast to the conical jet (A = 1), because the 
prefactor (2A — 1)~^(1 — -Tf /rf) can be larger than unity due to the slow acceleration 
r (X r^. The space-time diagram in Fig. [5] is useful to understand the following. 

(A) The radiation-dominated jet can cause internal shocks even during the acceler- 
ation if the jet is more collimated than the conical shape, A < 1 in Eq. ([5]). 

(B) The shock can occur at a large radius since an accelerating shell asymptotically 
approaches a light path (see dashed lines in Fig. [2]) . The shock radius is large 
when the interval At gets close to the critical time Atcrit where the time-fraction 
vanishes Tlx = I - At /Atcrit = in Eqs. dZD and ([8]). 

(C) The slow shell can be repeatedly shocked by the successive rapid shells, which 
can be ejected with shorter intervals ~ ro/2c-r2 than the slow shells ~ rQ/2cr^. 
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In the weak confinement, the causahty is lost across the jet {OF > 1). Although the 
angular structure of F would arise depending on the initial condition, the acceleration 
is likely reduced because of denser streamhnes (smaller V) than in the conical case. If 
the confinement is strong initially (inside the star) and weak later (outside the star), 
the oblique shock does not necessarily accompany the weak confinement phase. 

§ Jei confinement by cocoon - The jet would be confined by the progenitor star (long 

GRB) or the surrounding dense matter or wind^ (short GRB).'3Sl'''33~Sni' Note that 

the merger of two neutron stars ejects baryons before the collapse to a black hole 

even with an ordinary magnetic field strength.'33 As the jet runs into the stellar 

envelope, the shocked jet and envelope flow sideways (if < Fj^^) and inflate a 

cocoon, whose pressure confines the jet. We can roughly estimate the index A in 

Eqs. dSI) and ([6]) by considering the pressure balance at 

(I) the jet head: The longitudinal balance of the ram pressure between the jet 

~ Lj/irO'^r'^c and the stellar matter ~ p*c^lifi determines the jet head velocity 
c/3;,~c(L,/7rcV*r2)V20-i. 

(II) the cocoon-envelope boundary: The cocoon expands with a speed c/3c by balanc- 
ing the ram pressure of the ambient matter ~ /7^,c^/3^ and the cocoon pressure 
~ Ljt/ (c/Sfit) {cj3ctY which is the total energy deposited by the jet divided by the 

volume of the cocoon, so that /3c ~ {Lj/c^p^.r'^Y/^l3ll^ ~ {Lj/c^p^.r'^)'^/^9~^l'^ 
with r ~ cPht. 

(Ill) the jet-cocoon boundary: The cocoon pressure ~ P*c^li1 confines the radiation- 
dominated (hot) jet with the transverse pressure ~ Lj /irO'^r'^cF'^ , leading to 

d ~ (L,/c3p,r2)l/6r-4/3. 

Substituting the envelope profile ~ r~", F (x la Eq. ([5]), 9 oc r^~^ in Eq. ([6|), 
and a constant luminosity Lj ~ const, we can equate the radial dependence to find 

A ~ (e oc r("-i°)/i^ inside the star ~ r"") (9) 

The ideal radiative envelope has an index n = 3 (i.e., A ~ 1/2), whereas more 
realistic presupernova models such as 16Ti^ have shallower slopes n ~ 2 (i.e., 
more confinement A ~ 3/7). Because of the strong confinement A < 1/2, the jet is 
basically dissipative before the jet breakout. Similar dissipation is observed in the 
simulations.'^''^'S3 The Lorentz factor of the jet at the breakout is around 

rb ^ (^-^^ ~ 10-30, (10) 
Vnsco/ 

where ~ lO^'' cm is the stellar radius (breakout radius) and nsco ~ 10^ cm is the 
central engine size. The variability timescale is about ~ r^f/cF"^, which is ~ nsco/c 
for A = 1/2, so that the central engine timescale may be marginally preserved. 

Even after the jet breakout, the jet would be confined to some extent. A plausible 
source of pressure is the fast cocoon that just escapes from the jet head sideways at 
the breakout. The fast cocoon is thereby fast (ranging from F ^ 1) up to the jet 
Lorentz factor at the breakout ~ /], in Eq. (|10|) . and confines the following jet with 
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Fig. 3. Radiation-dominated internal shocks (RDISs) can occur if the jet is confined in 6 oc r^~^ 
with A < 1, for example, by the fast cocoon pressure [Eq. because the jet acceleration 

slows down r oc [Eq. ([5])]. The RDISs accelerate photons and create e*. The nonthermal 
photons are generated from extended radii, where lower energy photons originate deeper with 
more scatterings [Eqs. (|15|l and ()16|l ] before emerging from the photosphere [Eq. H14[) ]. In 
this paper, we do not account for the low-energy spectral slope. 



a lag time 6t over a long distance ~ c6tr^ ~ 10^'^ {6t /I s)(r'b/30)^ cm. Also, the 
transverse pressme of the fast cocoon is initially comparable to that of the jet pj, 
and then declines slower pc oc T"^ oc V'~'^^^ oc r~^/^ than pj ~ Lj/nO'^r'^cr'^ oc r~'^^ 
if the jet were conical (A = 1), because the fast cocoon becomes matter-dominatecj*^ 
with the coasting expansion V' ~ A-kt'^FA oc r^. For pc ~ Pj, the jet is confined, 

2 

A ~ -. (0 a r~^/^ outside the star) (11) 

The confinement gets weaker than inside the star (with index larger than the bound- 
ary A > 1/2) but still continues (in the weak regime 1/2 < A < 1). Thus, the 
dissipation proceeds via the internal shocks caused by the variability due to the cen- 
tral engine or the jet-envelope interaction. The initial opening angle would be the 
inverse of the Lorentz factor ~ -^f,"^ ~ 0.03-0.1 in Eq. (jlOp at the stellar surface 
To ~ r,, lO^'^ cm, and the product FO oc r^'^~^ grows afterward, as required for the 
afterglow observations of jet breaks. 

§e^ photosphere and high F with internal shocks - When the rapid jet catches up 
with the slow one, an internal shock forms, which dissipates the relative kinetic 
energy. The relative Lorentz factor at the collision is the same as the initial value, 

*•* The fast cocoon may be initially radiation-dominated and expand to _r Tp-ffc ~ after 
converting the radiation energy into the bulk motion, where 7p ~ 71, is the random Lorentz factor 
per proton at the shock, as in the coUisionless bulk acceleration.'^ Then, the Lorentz factor of the 
fast cocoon may range from _r ~ 1 to ~ . This effect could alter A in Eq. (|lip . 
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since F (x r"^ for both jets. For the radiation-dominated jet, the kinetic energy 
is almost carried by photons, which is transferred by Compton scatterings of e^, 
forming a radiation-mediated shockPS A smah fraction of photons are repeatedly 
upscattered across the shock, extending a power-law spectrum to high energies like 
the Fermi acceleration (Blandford-Payne mechanism).^- Even for a moderate 
optical depth, tt ~ a few, the Compton y-parameter y ~ '^t/3i2-^i2 ^^^^ exceed 
unity for the relativistic shocks. After the shock passage, the turbulent motion 
induced by, for example, Richtmyer-Meshkov instability, also provides the scattering 
centers, continuing the photon acceleration.^^ Then, an appreciable fraction of 
the relative kinetic energy can go to high-energy photons and create pairs. 

The pair creation can continue to a large radius since the collision radius 
rapidly grows as the jets expand owing to radiation pressure [see Eqs. d?]) and ([8]) 
and Fig. [2]. The jet acceleration holds as long as the pairs trap the radiation. 
The maximum Lorentz factor is limited by the condition that the system becomes 
transparent tt ~ n'^axct'^^^ ~ 1, where n'^ = f±L/A'Kr'^m(.c'^r'^ is the comoving 
number density of pairs, t'^^^ ~ r/cF is the comoving dynamical time, L is 
the isotropic total luminosity and f± is the energy fraction of pairs. Then, the 
coasting Lorentz factor is 



y 47rmeC^ro ) ^^^^^[sOj llQiOcmJ VlO^^ergs-ij ^^^^ 
at the photosphere of radius 

. (z±y'' . 1 . i„.3 (blY (^y" ( J-"- X' (14) 

"\FoJ \30 J VlOiOcm/ ergs-iy ^ ' 

where (and hereafter) we suppose a jet breaking out of the stellar surface, ro ~ ~ 
10^° cm and Tq ~ A ~ 10-30 in Eq. ([TOD, with A = 2/3 in Eq. ([nj in the last 
equalities. The pair energy fraction f± is about the energy fraction emitted above 
the pair creation threshold ~ F±me,(? ~ GeV, and /± ~ 1 is implied in Fermi bursts. 
The high Lorentz factor F± > 10^ is also consistent with the Fermi observations. 

^Nonthermal photospheric spectrum - The radiation-dominated jet emits almost all 
the energy from the photosphere, which is nonthermalized by internal shocks. The 
spectrum would have a thermal peak with a high-energy power-law tail. We illustrate 
the spectral shaping using a constant-luminosity model: equal-energy rapid jets are 
ejected after a slow jet at intervals of the causal time Atsh ~ '^o/2ci~'2 fo/'^cFf) 
(see Fig. [2]). The internal shocks can repeat in the dynamical time of the slow jet 
and create sufficient pairs to sustain the opacity. Thus, the pair density is 
mainly determined by the annihilation as n^. 1/ aTcAt'^^^ after a time At'^Yi 
following a shock. The interval of shocks is redshifted in the comoving frame as 
Ad'^Yi ~ -^^^sh ~ FrQ/2cF2 . Then, the typical optical depth of the jet evolves as 

^ _ n'arct'^ ~ ^ LI£L_ qc — oc r^^^A _ -1/3 ^5^ 

TT n±aTCTdyn Fro/2cFi F^ ' ^ ' 
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Note that tt oc n'^t'^^^ oc V'^^r/F oc r~^'^^^ oc r^^/^ (steeper) if were conserved. 

Under the photosphere > 1, thermal photons carry the dominant energy 
at a (nonrelativistic) comoving temperature /ifp^ak ^ w-gC^ with a number density 

'^peak ~ /±^"'±"^eC^//lt'peak- ^'^i'^ ^Ool bath, 6^ COols by AT'^/T'^ ~ -^^peak/^^eC^ 

in a single Compton scattering, and becomes thermal in a time less than the dy- 
namical time, ~ {n'^^^y^aTc)-^\T'J ATl\ ~ t'^^J^/Tr < t'^^^. The nonthermal 
photons with v' > fpg^k ^^^o 1°^^ energy 

in a single scattering by the cooled e^. However, the cooling time exceeds the dy- 
namical time, icooi7 ~ (n'_|_(TTc)~"'^|z^'/Z\j^'| > ^dyn, at hv' < m^.c?' Itt- Therefore, the 
spectrum remains nonthermal below the observed frequency ^ Tv' (~ const), 

y ^ oc r^^-i ~ r\ (16) 

hTT 

even if the photons are generated below the photosphere r < r± with > 1. At 
Tt > 1, the photons are trapped since the diffusion time is tdiff T^/n'_i_aTC ~ TTidyn- 
The internal shock extends the spectrum to a flat power law, injecting a fair 
fraction of the shock energy at the frequency in Eq. (|16p . In the constant-luminosity 
model, the injected shock energy is proportional to the number of shocks ^dyn/^^sh' 
so that the emergent spectrum at the photosphere is 

^.F^ oc oc oc z.(l-2A)/(3A-l) _ ^-0.3^ (^7) 

with Eqs. (jlSp and ()16p . Interestingly, this is consistent with the high-energy Band 
spectrum of GRBs uF^, oc ^-^-^^'^-^^^ The nonthermal photons are generated from 
radii extending over several orders, where the lower energy photons originate deeper 
below the photosphere with more scatterings in Eqs. (jlSp and ()16p . 

Our model naturally explains the observations that the nonthermal energy above 
i/peak is comparable to the thermal energy because the radiation itself nonthermalizes 
the radiation. In the constant-luminosity model, the slow jets (mainly for thermal 
energy) and rapid jets (mainly for nonthermal energy) have comparable energies if 
they are ejected for a similar duration of the causal time of the slow jets ~ rQ/2cFi 
(see Fig. ED . 

The rapid jets could become transparent via annihilation before colliding 
with the slow jet. The released photons expand conically and only a part of them 
~ {0/9o)'^ oc r"^-^^"^ could shock the slow jet. The loss of shocked energy softens 
the spectrum as uFi, oc r-'^^^^ x r^'^^^ oc zyi/(i^3A) ^ j^-i .^j^j^ -gqg^ g^^^ p!7|) . 

which is still consistent with the observations. The angular structure of A [below 
Eq. ([8])] could also cause diversity to the high-energy index. The released photons 
could be absorbed by the rear jets with different streamlines, or by the surrounding 
(optically thick) fast cocoon, possibly developing the hollow cone-structured jet with 
~ rj, ~ 10-30 for the shallow X-ray afterglow (see below). The decoupled from 
the rapid jet could be caught up with by other rear jets and create further e^. 

^Implications - 
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• GeV-TeV spectrum and delay - The photospheric spectrum may extend to fmax ~ 
r±r^2mec'^ ~ 10 GeV(r±/3000)(ri2/3)2,^ where we expect the creation cut- 
off and possibly the (blue-shifted and broadened) annihilation line|^''^ providing 
closure relation^^ for Fermi and CTA to verify the photosphere. We may also 
expect a spectral break at ~ r±ri2meCp due to the Klein-Nishina effect. 

The GeV onset delay observed in Fermi burstd^^^ may be explained by the 
leading jet that is not confined by the fast cocoon and, hence, cannot keep internal 
shocks and to a high Lorentz factor. In this picture, the delay time is about 
~ r±/cr^ ~ 1 s in Eqs. ([TO]) and ([H]), as observed. 

The (collisionless) internal shocks continue beyond the photosphere within 
outflows after the radiation decoupling, and could produce the observed extra GeV 
component, possibly up to TeV, by synchrotron or inverse Compton emission. 'l^''^ 
The luminosity can be appreciable up to ~ f±L. 

• Early afterglow - The prompt emission is radiatively efficient 1 — /± ^ 50%, and the 
remaining kinetic energy of powers the early afterglow. The energy fraction f± 
can be read from the prompt spectrum above the pair creation threshold ^ r±m(,(? ~ 
GeV. Thus, we predict a correlation between the early afterglow and the GeV- TeV 
prompt emission: a steep decay of X-ray after glox^^''^ accompanies a low GeV- 
TeV energy fraction. Further Fermi-Swift co-observations would be helpful. The 
shallow X-ray phas^'ll^ could be produced by the (matter-dominated) fast cocoon 
with r ~ 10-30 that receives energy via photons from e^-annihilated jets [below 
Eq. (|17p ]. The reverse shock emission would be soft because of the increased number 
of emitting leptons. 

• Cosmic Rays - The RDISs would produce few ultra high-energy cosmic rays and 
neutrinos,'^ consistent with the IceCube upper limit.^ The fast cocoon energized 
by photons from e^-annihilated jets might contribute to these emissions. 

• Central engine - The RDIS model requires less fine tuning of baryon load than 
previously thought, expanding the parameter space even to the original pure-leptonic 
model. '21' 'SI The baryon load could be prevented by the dipole field near a black hole, 
which may be rather necessary for the Blandford-Znajek effect to launch a relativistic 
jet despite nonaxisymmetric turbulence.'^' The ejected Poynting-dominated jet is 
strongly confined in the star or surrounding matter, possibly dissipating into the 
radiation-dominated jet via turbulence. 

• Unresolved problems - We leave the low-energy spectral indexP''21' uFi, oc v and 
the Amati and Yonetoku spectral relation^l2i>[I2i for future discussion. Since FT' = 
const, the peak energy reflects the initial fireball temperature. 
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